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The method of two-point quasiclassical Green's function is reviewed and its 
applicability for description of multiple reflections/transmissions in layered 
structures is discussed. The Green's function of a sandwich built of super- 
conducting layers with a semi-transparent interface is found with the help 
of recently suggested quasiclassical method [A. Shelankov, and M. Ozana, 
Phys. Rev. B 61, 7077 (2000)], as well as exactly, from the Gor'kov equa- 
tion. By the comparison of the results of the two approaches, the validity of 
the quasiclassical method for the description of real (non-integrable) systems 
is confirmed. 



1. INTRODUCTION 

One of the most efficient tools in the superconductivity theory is the 
method of the quasiclassical Green's functioroBBo. In recent decades, the 
quasiclassical technique has been successfully applied to a broad variety 
problems. Efficiency of the quasiclassical theory is due to the fact that 
it allows one to get description directly on the large spatial and temporal 
scales relevant for superconductivity, completely eliminating from the the- 
ory short "quantum" scales given by the Fermi wave length. Ignoring fine 
features on the quantum scale, the theory gives a closed set of equations for 
smoothly varying envelops. 

The quasiclassical approximation is applicable in superconductors be- 
cause the order parameter A varies on the coherence length £ ~ ■up/A, 
which usually much exceeds the Fermi wave length Xf ~ h/pF, vf and pf 
being the Fermi velocity and momentum, respectively. As first discussed by 
Andreevi, the slow spatial modulations of the order parameter cannot sig- 
nificantly change the momentum. Then, the two-component wave function 
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ip can be presented as if) = efi PF ' r (") where u and v are functions which 
slowly vary along a straight line, a classical trajectory corresponding to the 
momentum ppu. After the factorization, the Fermi wave length does not 
enter the theory anymore. 

The quasiclassical condition is violated if the potential energy varies 
on a scale of order of Xf e.g. when disorder (impurities) or interfaces are 
present. Disorder does not pose any difficulty as long as one is interested only 
in the impurity averaged properties: On average, disorder enters the equa- 
tions via the self-energy, acting like a smooth effective potential. Scattering 
by an isolated interface also does not create conceptional difficulties to the 
quasiclassical theory: one can include the reflection and transmission events 
by an appropriate boundary condition for the Andreev wave function □ or 
the Green's function Q. It is worthwhile to notice at this point that the 
quasiclassical scheme does not give correct results for quantities local in the 
p — r phase space. For instance, a particle with the momentum Ufx reflected 
from a wall creates the density, |^| 2 = \e tkpxX + re~ tkFxX \ 2 , which oscillates 
on the Fermi wave length due to the interference term 2Re re~ 2lkpX . How- 
ever, the interference pattern, which is very sensitive to the value of kp x , 
i.e. the direction of incidence, disappears after integration in a small region 
of the incidence angle. Although the quasiclassical technique ignores the in- 
terference of the incident and reflected wave, it gives correct predictions for 
coarse-grained averages: Observables like e.g. the electric current density, 
are given by the integral over the momentum space, and the low resolution 
description suffices. 

In many interesting cases such as domains in the high-T c oxides or arti- 
ficial mesoscopic structures, a few partially reflecting interfaces are present. 
In a multi-layer geometry, the quasiclassical method meets severe difficulties 
as noted in Ref. ||. The purpose of this paper is to identify the physics behind 
the difficulties and check the validity of the recently suggested quasiclassi- 
cal scheme by a numeric comparison with exact results for a double-layer 
structure. 

We attribute the difficulties in the theory to the existence of interfering 
paths formed by sequential reflections in a multi-layer structure. Consider, 
for instance, a sandwich built of two superconductors, left and right, of finite 
thicknesses, a"' r ) with a partially transparent interface, Fig. |l](a). We assume 
that the outer walls and the interface are ideal so that all the reflections 
and transmissions are specular. A fragment of a classical trajectory for a 
symmetric sandwich is shown in Fig. |l](b). 

Various processes of the Andreev reflection of an electron-like excitation 
coming on the knot A are shown in Fig. |2[ In the present context, the most 
interesting is the process shown in Fig. ^(c) where the particle makes a loop. 
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Fig. 1. (a) Sandwich build of two superconductor layers with the order 
parameters A^ ,r ) and thicknesses a^ l ' r ^ separated by a partially transparent 
interface, (b) The semiclassical wave which corresponds to a particle hitting 
the interface at the point (knot) A, is split by the interface into the left and 
right beams. If the left and right layers have the same thickness, the beams 
meet each other again at the knot B after the reflection by the outer walls. 



The amplitude of the process, unlike the processes in Fig. |2](a),(b) contains 
the phase shift exp[ip F (£ (r) - £W)] where and C® are the lengths of 
the left and right arms, respectively. As a result, the total amplitude given 
by the sum of various contributions, non-trivially depends on the value pp. 
One comes to the conclusion that the Andreev factorisation, rp = es Pf " r ("), 
which is the essence of the quasiclassical method, becomes impossible. 

It is clear that the loop contribution is very sensitive to geometry and 
surface roughness. Indeed, (independent) variations of £('' r ) on the scale of 
Xf destroy the interference. Besides, loops exist only if the trajectories which 
are split by knot A in Fig.|| meet again on knot B. This may happen only 
in an ideal structure; any imperfection on a scale ^ Xf would not allow this 
"accidental" event to happen. (In other words, any real physical system is 
non-integrable.) In a real system, multiple reflections/transmissions lead to 
the formation of single-connected tree-like trajectories as has been discussed 
in Ref. ||. The quasiclassical technique is expected to be applicable if one 
solves the equations on the tree rather than on straight line trajectories. 

Below, we compare exact and quasiclassical theories. For definiteness, 
we consider the trajectory resolved density of states (DOS) obtained by 
both approaches. It is clear from the very beginning, that agreement may 
be present only on a low resolution level: strong microscopic variations of 
the exact DOS in the p — r space must be smoothen by integration in a small 
region of the phase space. Even then, disagreement is expected. The point 
is that the exact calculations are manageable only for an ideal geometry 
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(a) (b) (c) 

Fig. 2. (a) and (b) Examples of processes included into quasiclassical theory. 
An electron is Andreev reflected back as a hole and follows the original 
trajectory. The loops are not formed in these cases, (c) The simplest example 
of a loop-forming process. An electron is scattered in the right arm and then 
Andreev reflected as a hole in the left arm. Such a process is not included in 
the quasiclassical approach which considers only tree-like (simply connected) 
trajectories. 



whereas the quasiclassical theory is expected to be valid only when some 
roughness is present. In the ideal geometry, there are always loops which 
survive averaging in the phase space, and their contribution may spoil the 
agreement. 

The paper is organised as follows. In Sec.^] we review the quasiclas- 
sical approach to the multi-layer systemsB. We introduce the trajectory 
Green's function and its factorization. The Green's function in the double 
layer system is found using both quasiclassical theory (Sec.^) as well as 
Gor'kov equations (Sec^j). In Sec.[|] the exact approach is compared with 
the quasiclassical theory. The applicability and the limits of quasiclassics 
are discussed. The results are summarised in Sec.^j. 



2. QUASICLASSICAL TECHNIQUE 

After EilenbergerB, the quasiclassical Green's (retarded) function <?^(t*), 
which is a 2 x 2 matrix in a spin-singlet Cooper pair superconductor, is 
introduced as 

§*(r) = l -[ dZG*(r) , (1) 
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where G p (r) is the Green's function G(ri,r 2 ) in the Wigner representation 
G R (r) = jG R (r + ^r-^)e^. (2) 

The Green's function for a given direction of the momentum n obeys 
a differential equation which couples the spatial points along a straight line 
parallel to n (or, more precisely parallel to the Fermi velocity v p correspond- 
ing the selected point of the Fermi surface in case of anisotropic spectrum). 
The line is the classical trajectory of the electron (hole). 

Alternatively, to Eq. (|]), one may use the formulation of the qua- 
siclassical technique in terms of the 2-point Green's function on classical 
trajectoriest3'cl. The 2-point Green's function obeys the following equations 

iv— + H R n (ri)J g R (xi,x 2 \n,R) = iv8(xi - x 2 ), n = R + xin , (3) 
g R (x 1 ,x 2 \n,R)(-iv-^+ H R n (r 2 )\ = iv8{x x - x 2 ), r 2 = R + x 2 n , (4) 

where the derivative in Eq. (||) operates backwards. The 2x2 traceless 
matrix H R n reads, 

frR 1R _ yR 

e,n ll e,n ^e^n i 

h? n =( £ ~ A r Ps ^ An ) , v = vn, (5) 
£ ' n y -A* n -e + v • p s J 

where is the order parameter (which may dependent on the direction n), 
and p s = — ^ A, A being the vector potential, and Y< R is built of the impurity 
self-energy and the part of the electron-phonon self-energy not included to 
the self-consistent field A. (Below, we omit R, n and e for brevity and use 
the notation g R (x\,x 2 )-) 

The propagator g R tends to zero at \x\ — x 2 \ — * oo, and g R is an analytic 
function of e in the upper half plane for any 21,2 including \x\ — x 2 \ = 00. 

Although the observables can be expressed via the quasiclassical 1-point 
Green's function (xi = x 2 ), the 2-point Green's function turns out to be a 
useful intermediate object. It gives a full physical description of the system 
in the approximation where part of the orbital degree of freedom is treated 
classically (no quantum broadening in the plane _L n), with a complete 
quantum treatment of the electron-hole degree of freedom. 

It is important that the construction based on the notion of smooth 
classical trajectories remains valid in the presence of disorder (or phonons), 
in the standard approximation when the scattering is included on the average 
via the self-energy (provided "ppl 3> 1, / being the mean free path). 
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2.1. Constructing 2-point Green's function 

The two-point matrix Green's function on the trajectory (n, R), is con- 
veniently built from column "wave functions" (ft, (ft = (") which satisfy the 
equation 

\v^- x + H R {x)\ <ft = 0, (6) 

here H R (x) stands for H R n (r) at the trajectory point r = xn + R. 
Denote (ft the row built from a column (ft by the following rule: 
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By virtue of the identity (H R ^j = ~ T yH R T y , the row (ft{x) built from 
a solution to Eq. (||), satisfies the conjugated equation 

iv£- + H R (x))=o. (7) 



dx 

It follows from Eq. (|6|) combined with (|7]), that 

^ (Mb) = . (8) 

This relation is valid for any pair of solutions <ft a and cfti. 

The Green's function is built of the regular solutions to Eq. @, i.e. 
solutions satisfying the following boundary conditions 

(ft + (x) -> , x +oo , , Q . 
— > , x —>■ -co . 

Denote </>±^ the normalized solutions, 

$ iN \x)<f> i +\x) = l. (10) 

The normalization is possible because the l.h.s. is a (finite) constant as it is 
seen from Eq. (||). 

The Green's function reads 

Indeed, it satisfies Eq. (j||) and Eq. (||) at xi 7^ X2, and is regular at \x%— a?2 1 — > 
00 by virtue of Eq. (|9|). Besides, the normalization in Eq. (10) ensures that 
the discontinuity at x± = X2, 

g R ( x + 0,x)-g R (x-Q,x) = i, (12) 

has the value required by the 5- function source in Eqs.(||), and (H). 
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2.2. 1-point Green's function 

To find observables like the electric current or charge density, one needs 
the Green's functions with coinciding spatial arguments, i.e. the 1-point 
Green's function. 

The 1-point Green's functions, g±(x) = g R (x ± 0, x), can be expressed 
via the normalized waves (see Eq. (]il"l)) 

= ^\x)^\x) , g R {x) = ^ N \x)4>f\x) . (13) 
This expression can be identically written in the form, 

9+^ = T ( U ( ^ +(xH-(x) , g_(x) = - 1 > _(x)0 + (x) , 
(f)-{x)<j) + {x) </>_(x)0 + (x) 

(14) 

where the normalization of the wave functions cft± is arbitrary. 
Note the projecting properties: 

gM = ±g R , g±g* = o , Sp g£ = ±1 . (15) 

Tagging electron- and hole-like excitations in accordance with the direc- 
tion of their propagation (±x directions) and considering examples, e.g. the 
normal state, one concludes that g R can be identified as the (quasi) electron 
part of the Green's function, and g R is the (quasi)hole one (and vice versa 
for g£). 

Denoting 

a = - , 6E^t, (16) 
V- u + 

where u± and v± are the components of <j)±, 



4>±(x) 



u±(x) 
v±{x) 



Eq. (14) becomes 



(17) 

As shown in Ref. |l0], the 1-point Green's function of the quasiclassical 
theory ("£— integrated"), g R , is given by 

g R = g R + g R . (18) 



i.e. 



g R = ^ ( . N) + 4>i N ¥ + N) ■ (19) 
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In terms of g R , 



^(9 R ±l) , (20) 
and the relations in Eq. (|l5|) lead to the well-known normalization condition 



& 2 



(4«) 



2 

= 1 



and 

Spg R = . 

Combining Eqs.(|l8|) and (|l7|), one gets 





(21) 

This parameterisation of the Green's function has been recently sug- 
gested by Schopohl and Makiliil (see, also Ref. |l^). The present derivation 
leads quite naturally to this decomposition, and clearly shows the physics 
behind it. Seeing that a and b may be interpreted as the "local" amplitudes 
of the Andreev reflection for electron and hole (see below), we call them the 
Andreev amplitudes. The amplitudes a and b obey the nonlinear Riccati 
equationEUti From the Bogoliubov - de Gennes equation for tp =J^), the 
Riccati equation for the ratio u/v has been derived by Nagato et ai±3. 



2.3. Knot matching conditions 

In the quasiclassical picture, particles move on smooth trajectories, usu- 
ally, straight lines characterised by the direction of velocity n (and the ini- 
tial position R). On interfaces, where potentials change on atomic scales, 
the quasiclassical condition is violated, and the quasiclassical wave function 
spreads from the original trajectory to those coupled by quantum scattering. 
Following Ref. |9], we call a "knot" the point where classical trajectories are 
tied together. 

In a general case, the knot ties together N in- and N out-trajectories. 
The in-trajectories (or channels) are those which have the the Fermi velocity 
directed towards the knot; the out-trajectory is characterised by the velocity 
directed from the knot. The in- and out-trajectories are somehow numbered, 
I = 1, . . . , N; we mark the outgoing trajectories with / so that k' stands for 
the k-th outgoing channels. 

On the quasiclassical scale ~ vf/A, the knot is point-like, and one 
can define the knot value of the trajectory wave function. Denote ipi the 
2-component wave function on the i-th in-coming trajectory, i = 1, . . . , N 
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at the point where it enters the knot, and analogously the knot value on 
the k— th outgoing trajectory. 

The knot matching conditions suggested in Ref. || read 

N N 

= S k'i^i ' VV = ^ S k'i^i > ( 22 ) 
i=l i=l 

where S^i are the elements of the unitary scattering matrix. In the spirit 
of the quasiclassical theory, S^n is the normal metal property taken at the 
Fermi surface; it is an electron-hole scalar. This relation generalises the 
matching conditions of Ref. |6] to the many channels case. 

Eq. (H) together with the matching conditions in Eq. ( p2|) allows one to 
find the 2-component amplitudes on trajectories with knots, and, therefore, 
the Green's functions. 



2.4. Matching Green's functions 

First consider an isolated knot mixing semi- infinite trajectories (with 
no more knots on them). With the origin chosen at the knot, the trajectory 
coordinate x n extends from — oo to on the ra-th incoming trajectory, and 
< xy < oo on the ^'-outgoing one. As before, the requirement, 

0_ OT (-oo) = , (j) +! k'(oo) = ,m, k = 1, . . . ,N , (23) 

uniquely defines the solutions 4>-^ n (x n ) and (p +i k'( x k') ( U P to a normalization 
factor). The knot values of the regular solutions are conveniently written as 

0_, m (a? m = 0) = , 0+^(^=0) = ^ ,m,k = l,...,N; (24) 

the parameters a m or are "bulk" properties independent on the knot. 
The problem is to find the knot values 



<f>+,i(xi = 0) = , 0_, n ,(x„, =0) = " , , l,n = l,...,N, (25) 




needed to evaluate 4>+^{xi < 0) and (f>- in >(x n i > 0) and, therefore, the 
Green's functions on the trajectories tied by the knot. On each of the trajec- 
tories, the Green's function in the immediate vicinity to the knot is expressed 
via the corresponding pair a; and 6; or a n > and b n i by the relation in Eq. (^l|). 

As has been shown in Ref. |9|, the boundary condition for the Andreev 
amplitudes a n i and b\ can be conveniently formulated in terms of the deter- 
minant 

V({a},{b}) = det l-Sa&b , (26) 
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built of the S-matrix and the diagonal matrices a = diag(ai, 02, . . . , ajv) and 
b = diag(&i/, b 2 ', . . . , 6jv')i a m and by entering Eq. ( |2~4| ) are bulk parameters 
insensitive to the presence of the knot. 

An important property of T> is that it is a linear function of any of a's 
and fo's, so that identically 

V = V$ ) +a l vP , V = V^' } + b n ,vf r) , n',l = l,2,...,N , (27) 



where 

V^=V\ ai=0 , V? = -^-V ,V^=V\ bnl=0 , V^ = ^L-V , (28) 



Also, T> is a sum each terms of which is a product of an equal number of a's 
and 6's. From this one concludes that 

YP$ = YPW. (29) 
I n 

These identities are useful for transforming various expressions. 

Rephrasing the procedure described in Ref. [|, the knot values of the 
Andreev amplitudes in Eq. (|2~5|) can be found by the following formulae: 



a «' = -^feo > h = ~wy (30) 



For the case N = 2, equivalent relations have been derived in Ref. 15 



Now, one is able to build the Green's functions Eq. fll7| ) on the l-th in- 
and n-th out trajectories. 



K' , -1 



(31) 

With the help of the identities in Eq. (|27|) and Eq. (|20|) , one can present 



Eq. (31) in the form, 



= / 2 ^ " 1 " 2 ^ a * ] Afl(n') = / 2 7n / - 1 ^(1 - 7 „0 > 

9 " ^ -|(l-7/) "(27/ -1) J ' 5 "^2 7n ,6 n , _(2 7n ,-l) y 

(32) 

where 

11 = , ln' = -^~. (33) 



It follows from Eq. (j29j) that 
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This identity is directly related to the current conservation. 

Summarising, the Green's functions on trajectories linked by a knot are 
calculated as follows. First, one solves Eq. @ with the boundary condition in 
Eq. ( |23| ) on each of the trajectories and calculates "bulk" functions <p- trn (x < 
0) and <f) +i k'(x > 0). By this, one finds a m and by in Eq. (^) . For a 
knot characterised by a scattering matrix S, one is then able to find the 
determinant T> Eq. (^). The next step is to calculate the knot values of 6's 
on the incoming trajectories and a's on the outgoing ones using Eq. (j30|). 
Then, one solves Eq. (||) for </> +j / (x < 0) on the incoming trajectories and 
4>- t n'(x > 0) on the outgoing ones with the boundary condition <^ +j /(0) = Q) 
and (f)- >n r(0) = yf'), respectively. The 1-point Green's function is then built 
from (fi± using the representation in Eq. (p]). The knot values of the Green's 
functions can be also found from Eqs.(|31|)7or (|32|). 

This scheme is also applicable when the trajectories connected by the 
knot under consideration may enter other knots, in a tree- like trajectory (see 
Fig. ||). As a matter of principle, one assumes that the system under consid- 




Fig. 3. An example of a tree-like trajectory. Pieces of the straight lines show 
the trajectories before or after they enter a knot (filled circles), i.e. before 
or after a collision with an interface. There is only one path connecting any 
two points on the tree so that the tree is effectively 1-dimensional. 

eration is finite, and it is surrounded by a "clean" material where trajectories 
are infinite lines without knots. Then, one solves the problem for the knots 
on the boundary and moves inwards towards the knot of interest. In the 
one-dimensional topology of the tree with only one path connecting any two 
knots, the procedure is unique. 
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3. SANDWICH: QUASICLASSICAL GREEN'S FUNCTION 

In this section we calculate the quasiclassical Green's function in a dou- 
ble layer system depicted on Fig. [j]. The real space classical trajectory of a 
particle in a two layer system is formed by multiple reflections on the outer 
surface and the interface between the layers (see Fig. |].). 

The Green's function can found using Eq. (|32|). In this simple case with 
only two in- and out-trajectories, determinant T> reads 

V = R{1- ai b v )(l - a 2 b 2 ,) + T(l - ai6 2 /)(l - a 2 b v ) , (35) 

where R and T are the reflection and transmission coefficients. 

The problem is to find the amplitudes a m on the incoming legs of a 
certain knot and by on the outgoing legs. Since none of the legs lead to 
infinity the situation is slightly different than in the previous chapter and 
the procedure has to be modified. 




Fig. 4. The trajectories of a particle in a two layer system. 



The argument used to find the equations for the amplitudes a m and by 
is the invariance with respect to translation parallel to the interface. Suppose 
we know a transfer matrix which relates the wave function </>_ in the legs 
say " 2" and " 4" . Because of the translational invariance these two functions 
can differ only by a prefactor. In another words (j) ls eigenfunction of the 
transfer matrix and the ratio of the components of (ft must be the same in 
both corresponding legs. Thus the amplitudes a in legs "2" and "4" are equal 
a 2 = 04. One can argue in the same way to show that a± = 0,3, by = by and 
so on. 
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We can proceed in the following manner. First we find the transfer 
matrix as a product of knot transfer matrix and propagator across the layer. 
Then we relate the amplitudes in corresponding legs. We get a set of 4 
equations for 4 unknowns ai, 02,61' and by (see Fig. ||. for notations). 

The knot transfer matrix connecting the wave function in the incoming 
and outgoing legs on a knot can be written as follows 

*"«u?)(«!)u?r- (36 » 

The transfer matrix!! for the channels on the right side is a function of the 
amplitudes only on the left side. 

The propagator across the layer is found by integration of Eq. (|6|) . For 
homogeneous superconductor the expression reads: 

/ M \ / i^- 

ut \ 1 a o H e F 

W(a;) oc (r) " w 

V 6 o W( e~ l — 

Here, aj r) = A< r ) /(e + f W) and = A( r )*/(e + ^) are the bulk values of 

a and b in the right layer, and = \J e 2 — |A( r )| 2 (Im £( r ) > 0). 

Once the transfer matrices for both sides are known, one can write down 
the 4 equations for the four unknown amplitudes axp an d £>i'2'- Solving 
the equations and using the expressions relating a's and b's in the in/out 
going channels from the previous chapter one finds the amplitudes in every 
leg. Thus the Green's function can be constructed and the density of states 
calculated. 

When the numerical calculations are carried out, it is easier to use an 
iterative procedure instead of solving the set of four nonlinear equations. 
One starts with the bulk values of a and b in, say, right side and calculates 
the corresponding amplitudes on the left. These values are then used to 
calculate back a and b on the right side. This process is continued until 
self-consistency is reached. 

4. SANDWICH: EXACT GREEN'S FUNCTION 

In this section we calculate the exact Green's function for a sandwich. 

Consider a sandwich formed by two homogeneous layers with thicknesses 
a G> r ) anc [ order parameters A('' r ); here and below, the superscript (I) and 
(r) label quantities which refer to the left and right layer, respectively, see 
Fig. |(a). 
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The left (right) layer occupies the region — o'^ <x<0(0<x< a^). 
The interface is characterised by the scattering S'-matrix, 

s =(l r t )- < 38 > 

with the transmission and reflection amplitudes t and r for the case of 5- 
function interface potential. The outer walls are impenetrable. The goal is 
to find the exact Green's function G(r,r') from the Gor'kov equation. 

Due to the translational invariance in the direction parallel to the layers, 
the problem is effectively one-dimensional. For a given value of the parallel 
momentum py, the Green's function, G = G £tP ,,(x,x'), depends on the two 
coordinates x and x' as well as the energy e. The Gor'kov equation reads: 

(e - H)G(x, x') = I6(x - x') ; H = ( J, A ^ . (39) 

For the parabolic spectrum, £ = p 2 /2m—p 2 Fx /2m, where pp x is ^-projection 
of the Fermi momentum p 2 Fx = pp — p?,. The Green's function is continuous 

at x = x'\ G(x + ,x) = G(x~,x). Its derivatives suffer a jump generated by 
the delta- function on the r.h.s. of Eq. (|39"|): 



x=x'+0 

p x G(x,x') =2m/ih. (40) 

x=x'— 

The boundary conditions at x = corresponding to the semi-transparent 
interface are conveniently written as 

P a G(0 + ,x')= Maa'P a 'G(Q-,x') a = ±, (41) 



vhere the projectors P ± , 



2 

and the transfer matrix A4, 



o'=± 



P ± = l -(l±p x /p Fx ) , (42) 



l/t* r/t 



M = [ r*/t* l/t ) ■ (43) 

It is convenient to express the Green's function in terms of two-component 
functions $> V)(T (x) and <& u ^{x') defined as follows: 

$„a{x) = fae^"* ; $ ua (x') = 0e- l(7p » x ' . (44) 
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Here a = ±, p v = pf x + v^jm is the electron (y = +) or hole (y 
momentum, and <j> v is the solutions to matrix equation 

A* -tf ) ^ = ^ ' 



(45) 



and (j) u is the adjoint row- vector defined in such a way that (j) u (j)i_ l = 5 Ufl . 
The Green's function can be then written as: 



ihvF x G(x, x')a z 



E $ v , tT (x)(v,(r\Q > \i/(T , )$ l/ > jtT >(x') forx>x' 
E ^(xX^ctIQ^z/ct'^vKx') for x < x' 



(46) 

The 4x4 matrices Q <>> are constants. It has different form for x > x' and 



x < x'. 



Inserting the Green's function into Eqs.( |3 9^400 and using the match- 
ing condition at the interface in Eq. ( [4T| ) one gets (u, a\Q < ' > \u'a'}. After 
straightforward but lengthy calculations, we have found the following rather 
simple result: 



(47) 



[ TzQll TzQ\2 \ _ J_ [ W + ~ A + A_W + W- A-W + (W- - 1) 

y zQ 21 T zQ 22 ) D \ A + w-(w + -l) W^-A + A^W-W + 

D = 1 - A_A + Ti(w + w^) . (48) 

Here, w± is a 2 x 2 matrix, w± = W±/TiW±, 

^-(^W' + MV^P + «"•*'•<'>) , (48) 



and A± are (scalar) coefficients, 



which have the physical meaning of the Andreev reflection amplitudes for a 
transparent interface. 

The Gor'kov Green's function for a sandwich is known from the liter- 
ature 0. Eq. (|47|) agrees with the previous work, giving simple and concise 
form for the Green's function. 
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4.1. Coarse-grained Green's function 

As have been already mentioned in the Introduction, the quasiclassical 
theory does not attempt to give any good description on the truly microscopi- 
cally local level. Instead, it supplies knowledge about coarse-grained observ- 
ables. Accordingly, one should derive the coarse-grained Green's function 
corresponding to the exact theory before comparing it with the quasiclassi- 
cal counterpart. 
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Fig. 5. Discrete energy levels in then sandwich Fig. |l|(a) corresponding to 
eigenstates with a definite momentum pn. The curves show how the position 
of the level e changes under small variations of the angle 0, cos G = P\\/pf- 
The sandwich with a"' = up/A^ is thick on the scale of the Fermi wave 
length: aft' = W00/pf] the transparency of the interface T = 0.9. Other 
parameters of the sandwich: = 3a^ l \ = — 2A^. 



The space averaging of the partial (i.e. momentum resolved) density of 
states is rather simple. For a given pu, the density of states is proportional to 
the imaginary part G p ., (x, x' — > x). It is seen from Eq. (44) that the elements 
of the Q-matrix which are off-diagonal in a, create rapidly oscillating terms. 
Therefore, space averaging amounts to ignoring the off-diagonal terms. 

Even reduced to a low spatial resolution, the Green's function remains 
to be a fast function of the momentum direction. For instance, the density 
of states is non-zero only at those discrete values of pu where the energy 
variable e equals to the energy of the bound state for the finite motion 
in the x-direction. Fig. || shows the fast periodic motion of the energy level 
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positions under a small variation of pit . An integration within a small interval 
of the angles is required to transform the density of states from a sequence 
of the 5-function spikes into a smooth function of the momentum direction. 

3.5, r— i 1 , , , , , 1 
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Fig. 6. The energies of the bound states as a function of angle (full 
line) superimposed on the quasiclassical angular resolved density of states 
(dashed line). The parameters of the double layer are A^ r ^ = — 2A^',a^' = 
v F /^ l \a^ = 3u F /A«,T = 0.9,G = vr/4. 

In the next section we compare the coarse-grained density of states 
extracted from the exact Green's function with the quasiclassical trajectory 
resolved DOS. 



5. DOS: QUASICLASSICAL VS. EXACT 

In the most simple case where the order parameter is the same in the 
left and right layers, the both theories reproduce the BCS density of states 
for any transparency of the interface (as a consequence of the Anderson the- 
orem). Nontrivial comparison requires an inhomogeneous order parameter. 

First, we analyze a generic "sample", parameters of which are in no 
special relation to each other. In Fig. |(| the energy of the bound states of 
the exact theory are superimposed on the quasiclassical DOS. 

One sees that indeed the periodic motion of the levels occurs in the 
allowed bands predicted by the quasiclassics. From Fig. 0, which shows the 
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DOS as a function of energy, one may conclude that the quasiclassical theory 
gives correct energy dependence, excepting perhaps some very fine details. 
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Fig. 7. The angular resolved density of states calculated using both exact 
(dashed line) and quasiclassical (full line) methods. The parameters of the 
double layer are A W = -2AW jfl W = v F /A®, aM = 3v F /A®, T = 0.9, 6 = 
vr/4. 



In Fig. ||, we show the DOS for a specially selected "sample" for which 
non-classical effects are expected to be most pronounced. The sample is left- 
right symmetric except that A® = — A^ r ^. Because of the symmetry, the 
contribution of the loop in Fig. |2|(c) survives the coarse-grained angular aver- 
aging. Besides, the knot scattering matrix S of the "sample" is intentionally 
chosen as i.e. 

• * 2 -(?o). 

In this case due to the quantum interference of reflected and transmitted, 
the interface is fully transparent if the particle hits it twice (provided the 
sandwich is symmetric) . The quasiclassical theory misses this phase sensitive 
effect since only the probabilities enter the theory. 

As clearly seen from Fig. ||, the two approaches show distinct result 
for this "sample". We attribute the disagreement to the loop contribution. 
This point of view is supported by Fig. |9| where DOS is shown for the sample 
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Fig. 8. The quasiclassical (full line) and exact (dashed line) angular resolved 
density of states for V NOT case. The parameters of the double layer are 
A (r) = _A(0 )O (0 = 4v F /A ( - l \a^ =4v F /AV\T = 2/3,9 = vr/4. 

differing from the previous one only in a slight violation of the left-right sym- 
metry. For this sample, the disagreement is much less pronounced. This can 
be understood since in an asymmetric sample, the simple loop contribution 
averages to zero. 

6. CONCLUSIONS 

The main goal of the study has been to check the validity of the quasi- 
classical scheme suggested in Ref. || for the description of multiple reflections 
in layered superconducting systems. For this, we have compared the angu- 
lar (trajectory) resolved density of states obtained from the full Gor'kov 
equation and the quasiclassical theory. From the comparison, we conclude 
that the approach based on the notion of quasiclassical tree- like trajecto- 
ries El is in qualitative agreement with the exact theory. Some quantitative 
disagreement which comes as no surprise, seems to be under control. The 
point here is that the assumption about the tree-like character of classi- 
cal trajectories can be justified only for rough interfaces when the system 
becomes non-integrable, whatever small the disorder. However, the exact 
Green's function is calculated in the ideal geometry where the trajectories 
have loops (as in Fig. |l]). We ascribe the disagreement to the loop contri- 
bution. Choosing the geometry, we able to control to some extent the loop 
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Fig. 9. The angular resolved density of states for "almost" VNOT case. 
The parameters of the double layer are A( r ) = — AW, aft> = 4vf/ 'AW, = 
A.lvp / A^ l \T = 2/3, = 7r/4. The full line corresponds to the quasiclassical 
theory, the dashed line to the exact calculation. 



contribution. Indeed, the disagreement is most pronounced in the symmet- 
ric "sample" (see Fig. ||) where the contribution of the simplest loops shown 
in Fig. |](b) survives coarse-grained averaging. A small asymmetry, which 
destroys their contribution, considerably reduces the disagreement (see Fig. 
^). Note that loops always exist in the ideal specular geometry: even for an 
asymmetric sandwich, a® ^ loop is formed after every second colli- 

sion with the interface. These higher order loops are not destroyed by either 
the angular averaging or the averaging with respect to the layer thickness 



as in Ref. 13 and 14. In our opinion, the difficulties with the quasiclassical 



theory reported in Ref. [L3| and 14 originate in the contribution of the high 
order loops. 

Finally, our conclusion about the validity of the quasiclassical theory is 
as follows. The theory is applicable for the description of a typical real sam- 
ple where some roughness of the interfaces and surfaces is inevitably present. 
Since real typical samples are most probably non-integrable, the trajecto- 
ries are tree-like, and the scheme suggested in Ref. |9| is valid. However, the 
scheme does not include certain physics which may be of importance in some 
situations where the sample is manufactured to enhance the role of interfer- 
ing paths. Then, one should use the full version of the Gor'kov equations 
for the description of such a special system. 
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